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Abstract The mechanism of the roughening induced partial depinning of glid- 
ing dislocations from 3 He impurities is proposed as an alternative to the standard 
"boiling off. We give a strong argument that 3 He remains bound to dislocations 
even at large temperatures due to very long equilibration times. A scenario lead- 
ing to the similarity between elastic and superfluid responses of solid 4 He is also 
discussed. Its main ingredient is a strong suppression of the superfluidity along 
dislocation cores by dislocation kinks (D. Aleinikava, et. al., arXiv:0812.0983 1. 
These kinks, on one hand, determine the temperature and 3 He dependencies of 
the generalized shear modulus and, on the other, control the superfluid response. 
Several proposals for theoretical and experimental studies of solid 4 He are sug- 
gested. 
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1 Introduction 
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The discovery of the torsional oscillator (TO) anomaly by Kim & Chan in 2004 1 , 
which was originally ascribed to the non-classical moment of inertia (NCRI) 2 , 
has reignited a strong interest to the supersolidity— in solid 4 He. However, after 
observing the enigmatic similarities between the TO and the shear modulus by 
Day & Beamish 4 (the similar effect has been observed by Tsymbalenko 5 at sig- 
nificantly higher temperatures and frequencies), serious doubts have been raised 6 
about the supersolid interpretation of the TO anomaly. Very recently in Ref . 7 it has 
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been shown that the deformation of solid 4 He in the torsion rod is fully responsible 
for the anomaly reported by several groups. 

A completely different approach has been undertaken by Ray & Hallock— who 
have observed the direct superflow through solid 4 He at very small rate - about 
3g/year. Thus, it is very likely that such a signal is below the sensitivity of the 
TO 8 . Under these circumstances, it is natural to ask: is there any "room" left for 
observing the supersolidity in the TO-approach? 

While not providing a definite answer to this question, we show that the su- 
perfiuid (SF) response of the dislocation cores (found to be SF in the ab initio 
simulations*^) can mimic the temperature and 3 He dependencies of the gener- 
alized shear modulus 4 (see Eq. d22b ). The key to such a similarity is a possibility 
that dislocation kinks strongly suppress the superfluidity along the core 11 (cf. also 
a proposal by Balibar— ). 

Our paper is organized as follows: in Sec|2]the shear modulus anomaly— is 
discussed in terms of the dislocation roughening in the presence of the Peierls 
and 3 He pinning potentials. The inconsistency of the 3 He "boiling off' scenario is 
revealed, and it is suggested that 3 He remains bound to dislocations even at high 
T. The puzzle of the "missing dissipation" is addressed in Secf3] Then, in Sec|4] 
we will discuss the coarse grained model where the SF order parameter interacts 
with the thermal kinks, and will show that the condensate fraction can resemble 
the shear modulus T- and the 3 He- dependencies. We also discuss a mechanism 
linking the dissipation in the SF and in the gliding dislocation subsystems to the 
common bath of thermal kinks. Finally, in Secf5]we give a brief summary of the 
results and outline proposals for future studies. 



2 Generalized shear modulus 

The shear modulus anomaly of many materials is caused by liberation of glid- 
ing dislocations from the low temperature pinning by either impurities 13,4 or by 
the Peierls potentia l n i 14 . As shown in Ref. 15 , the reconciling of the observed T- 
dependencies of the shear modulus G(T) with the "boiling off-model of 3 He 
impurities from the dislocation cores requires introducing a wide distribution of 
3 He activation energies E a - as wide as the mean energy itself. The origin of such 
a wide distribution in a crystal with relatively low dislocation densities ( about 
w 10~ 7 — 10 -8 in the atomic units of the typical distance between 4 He atoms 
b « 3.5 — 3. 7 A) is very puzzling and raises more questions than gives answers. 

In this section we argue that the actual shape of G(T) in the presence of the 
Peierls potential and the 3 He pinning is determined by thermal roughening of glid- 
ing dislocations: thermal and quantum fluctuations of dislocation shape wash out 
both potentials. This process is the alternative to the "boiling off' mechanism. 



2.1 Fluctuative depinning of gliding dislocation 

In solid 4 He dislocations are strongly pinned by the Frank's forest cross-linking 
points 13 and, in addition, are weakly pinned by Peierls potential and 3 He impuri- 
ties. The model capturing both effects relies on the Granato-Liicke-type string— 
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description ' subjected to both the Peierls potential Up and the trapping poten- 
tial U, provided by 3 He atom. The corresponding action capturing quantum and 
thermal effects can be written in the imaginary time % as 

H = J^dxj o dr (J^ ((V x y) 2 + (V z y) 2 ) -ay- U P (y) -N(x)U t (y)\ , (1) 

U P (y) = -acos(27Cy), U,(y) = N(x)V(y), V(y) = -V Q ^ V (-(y/y ) 2 ), (2) 

where y(x, t) is the dislocation displacement from its equilibrium position in a 
gliding plane; L stands for the dislocation length between the cross-linking points 
(y(x = 0,z) = y(x = L,x)=0); j3 = 1/T (in atomic units h = 1 , K B = l)\K denotes 
the Luttinger parameter 14 ; a denotes external stress; a,Vo stand for the strength 
of the Peierls and the trapping potentials, respectively; N(x) gives the density of 
3 He impurities which can be located at any point < x < L along the line y = 0; 
the trapping potential is chosen as a Gaussian with some range yo ~ 1 . 

In Eqs. (11121) and below all lengths are measured in units of Burger's vector 
b, and the unit of time is chosen so that the speed of sound V s along the core is 
unity. Accordingly, the unit of energy (temperature) is given by T„ = hV s /Kgb = 1 
(which is ~ 10K in standard units). All our results below are represented in such 
dimensionless units. 

The parameters can, in principle, be extracted from the ab initio simulations. 
For example, Ref. 17 provides the estimate Vq « 0.8K for binding of 3 He to screw 
dislocation. In our simplified model we ignore the long-range part of the trapping 
potential, which should scale as V(y) ~ — l/\y\ for the case of edge dislocation 
(and as ~ — \/y 2 for the screw dislocation—). 

In Ref. 18 it has been suggested that 3 He atoms don't actually provide pinning 
for dislocations, and, instead, they induce a viscous drag force. Here we take a 
different approach: due to the extremely narrow band width J ~ 10" 4 Ki2, a 3 He 
atom-impuriton, once bound to a dislocation, can be easily localized in the lattice 
potential to become an immobile pinning center. 

Monte Carlo simulations of the partition function Z = 
/ Dyexp(— H) as well as any needed mean (...) = / Dy...exp(—H) /Z have been 
conducted in discretized space-time lattice using the gradient expansion method 
combined with the Worm Algorithm?- with the 3 He impurities considered as 
Boltzmann particles of fixed total number, with N(x) = 0, 1 depending on if 3 He 
atom is or is not present along a given space-time line (x, t), respectively. This 
model is similar to the one considered in Ref. 14 where Vb = 0. Here we ignore the 
long-range potential between kinks 14 . 

The response of the model (fl} can be related to the average coarse-grained 
shear modulus 14 . If y(x,t) is a typical dislocation displacement of an element of 
Frank's forest of typical sizes L,L y ,L z (along the respective axes X,y,Z), the av- 
erage strain is ta J y(x,t)dx/(LL y L z ) (in the chosen units). Since y a, and the 
total strain must include the elastic one ojG e i, with G e \ being bare shear mod- 
ulus, the resulting (inverse) shear modulus for static a becomes G _1 = G~} + 
A (y) j (LL y L z a), where A is some orientation averaging factor 16 ~ 1 . In the limit 
<7 — > this gives 

1 1 (jjA rP r L r L 

- = — + jL- dx dx dx>{y(x>,T)y(x,0)) (3) 

Cj U e i L Jo Jo Jo 
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Fig. 1 (Color online) Shear modulus G from Ref4 for 1 ppb of 3 He at / = 2000Hz (blue dots) 
and the MC results d 1 13b for various (linear) concentrations n of 3 He impurities. 
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Fig. 2 (Color online) The experimental data 4 are the same as in Fig.Q] The family of the MC 
curves from Fig. [T] has been collapsed to a single master curve by rescaling the T-axis and 
slightly adjusting G e i to allow all curves to have the common value at T = 0. The simulation 
parameters are chosen as K = 0.1, Vb = 0.3, y^ 1 = 4.2 and the Peierls potential amplitude a = 
0.01 (in units of T r ,,b) was adjusted to achieve best agreement with the experiment at n = 0. 

where = 1/ (L y L z ) w \/L? « 1 stands for density of dislocations in the units 
of b and the mean (...) is evaluated with respect to the action OJ at a = 0. 

The results of the simulations are presented in Figs. II I2| for the case of pure 
Peierls potential (« = 0, cf.— ) as well as for a set of finite linear densities n of 
3 He. For larger n the softening of G(T,n) occurs at larger T. Such tendency is 
clearly seen in FigfT] We have found that the curves for different n are almost 
self-similar to each other and can be collapsed to a single master curve, Fig.|2j by 
choosing a proper multiplicative factor 7^ : T — > T^T. The simulations have been 
conducted for fixed values of n, that is, 3 He atoms were not allowed to boil off 
from the dislocation. Despite that and the strong pinning at the dislocation ends, 
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Fig. 3 (Color online) The shear modulus (/ = 2000Hz) normalized by its maximum variation 
for two different 3 He concentrations 1 ppb (blue dots) and 0.3ppm (open red rhombi) 4 . The data 
for 0.3ppm is also shown for the rescaled temperature T — > Q.67T (filled rhombi). The Monte 
Carlo data (stars) for n = 1/40 of 3 He from Fig.[T]is shown at the same rescaled temperature. 
Black solid line: the attempt to fit the 0.3ppm data by the "boiling off model with the parameters 
chosen n d = 0.3 • 10~ 8 , E„ = 0.65K. 



thermal kinks induce fluctuations of the dislocation position y around y = so 
that the Peierls and the 3 He-pinning potentials are essentially washed out, and the 
modulus can reach its high-T value. 

The rescaling parameter reflects the value of the dislocation gap A = 
induced by the pinning potentials. We have found that A grows with n. This 
can be understood from the following consideration: at small n the effect of im- 
purities is inducing some effective average potential ~ — «Vo ex p(~ (y/yo) 2 ) ~ 
— nVb(l ~~ (y/yo) 2 ) in the action OJ. This corresponds to the formation of the ef- 
fective gap ~ A Fts y/nVo/yo- As n grows, the effect of the impurity-impurity repul- 
sion through the string excitations becomes important ( the pinning increases the 
zero-point energy of the dislocation phonons and kinks). Thus, the gap must show 

a crossover from A ~ A^ + ...n at small n to A ~ n at large n, where Aq ~ a is 
attributed to the Peierls potential. The fit by such quadratic function gives A (n) — 
v /A 2 {0) + din + 5 2 n 2 1 A (0) = 0.209 ± 0.005, ft = 1.17 ± 0.06, 8 2 = 5.3 ± 0.3. 

It is important that, while strongly fluctuating dislocation becomes practically 
free from the pinning potentials, it still provides a very strong binding potential Vq 
to the impurities on average if compared to the impuriton band-width /. Indeed, 
we estimate Vb as V(y), Eq.©, smeared by large fluctuations of the dislocation 
position y around y = as 



(4) 
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where the thermal fluctuations are taken in the high- 7 limit as (y 2 ) « TL/Td > y^ 
and To stands for Debye temperature. An estimate of Vq/J for the typical values 
(V ~ lK,r - IK, T D ~ 10K, y ~ 1, L ~ 1 j y/n~ d ~ 10 4 - 10 5 ) shows that %/J ~ 
10 2 . As discussed below in Sec l2.2l such large ratio implies very long equilibration 
time for 3 He to boil off into bulk. 



2.2 Fluctuative partial depinning versus boiling off of 3 He 

Figs. [T|2] represent the result of the fluctuative creep — quantum and thermal fluc- 
tuations of the dislocation shape due to kinks and phonons produce decoupling of 
the dislocation line from the Peierls and 3 He-pinning potentials. In other words, 
despite the presence of impurities at the equilibrium dislocation position (y) = 0, 
the spatially averaged pinning force acting on the dislocation is significantly di- 
minished by the shape fluctuations. The question is how allowing 3 He to boil off 
will modify the curves G(T,n). 

Let's, first, turn attention to the experimental situation. Fig. [3] shows Day & 
Beamish data for G(T) at two vastly different bulk concentrations of 3 He. These 
moduli can be collapsed on each other and fit by the Monte Carlo data obtained 
under the assumption that 3 He does not boil off. In contrast, the boiling off model 
with the single activation energy can fit neither of G(T) (cf.— ). In other words, 
the experiment indicates that 3 He remains bound to the dislocation even at high 
temperatures, which is consistent with the conjecture of very long thermalization 
time (as our estimates indicate — much longer than any reasonable experimental 
time) of the impuritons. 

The equilibrium concentration n(T) of 3 He on dislocations at density n c i can 
be found within the simplest thermodynamics consideration treating 3 He atoms as 
non-interacting Boltzmann particles^: 

n(T)= - X ' (5) 
n d + exp(-E a /T) 

where we assume the total (bulk) fraction Xt, of 3 He as Z3 < (in the chosen 
units) and E a stands for 3 He activation energy. This equation is also valid for 
X 3 > n d as long as T > E a /\ ln(X 3 - n d )\]. 

The boiling off starts at 7t, on f» E a /\ ln(n d ) \ «E a for the realistic densities of 
dislocations n d ~ 10~ 6 — 10 -9 . Given the typically accepted E a ~ 0.5 — 0.8K, this 
gives the boiling off onset temperature as 7t, on ~ 40 — 70mK. The pinning becomes 
irrelevant when the typical distance r ~ 1/n between the trapped atoms is larger 
than L rj 1/ s/nd, Ref. 4 . This determines the upper temperature 7i, up when all 3 He 
atoms have essentially evaporated as 7i, up ~ E a /\\n(X^/ ^/n2)\. For the standard 
values of X 3 ~ 10" 6 - 10" 9 << 1 this gives r bup w 70 - lOOmK implying that the 
actual softening range should be much narrower (< 30mK) (see the solid black 
line in Fig. [3) than it is observed experimentally—. 

The pure boiling off model corresponds to the action OJ without Peierls po- 
tential (a = 0). In the approximation of strong pinning (Vq — > °°) by the impurities 
placed (equidistantly a distance r = 1 /n(T), Eq.©, apart from each other) along 
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the dislocation of the length L « 1 / ^fn/i one finds (using the free string solution— 
between two pinning centers) the shear modulus (f3]) 

l J_ y = _i i_ 

G(r) ~ G e , + 1 + n 2 /n d ' 7 ~ G(~) G e7 ' 

This equation together with Eq.© have been used to obtain the solid fitting curve 
in Fig. [3] for the parameters typical for solid 4 He. As can be seen, this fit is not 
adequate^. 

Similarly, the attempt to include the boiling off mechanism into the Monte 
Carlo data by replacing n by the form (0 leads to the same failure if one chooses 
the realistic values of n,j,E a ,X^. Simply saying, if the equilibrium were reached 
fast on the experimental scale, there would be no impurities left on dislocations 
to affect the curve G(T) at T > 0.07 — 0.1K. Nevertheless, as discussed above, 
the fluctuative mechanism completely ignoring the evaporation of 3 He from dis- 
locations can describe well the temperature dependence of G(T) in a wide range 
of 3 He concentrations. In other words, it appears that the actual equilibrium be- 
tween binding to dislocations and being free in the bulk cannot be established for 
3 He atoms in any reasonable experimental time. Accordingly, the amount of 3 He 
trapped by a dislocation is predetermined during crystal growth. 

The following considerations may help to understand our conjecture: 3 He 
atoms in a solid 4 He exist as impuritons with a very narrow band-width J so that 
any repulsive or attractive potential characterized by energies larger than J causes 
a 3 He atom to be localized within a spatial region of a size on which the potential 
energy changes by not more than J, Ref. 19 . Furthermore, bulk phonons which are 
responsible for establishing the equilibrium can only change an impuriton energy 
by J which is much smaller than Vq, Eq.®- In other words, there is no chance 
for an impuriton to be freed into the bulk by just absorbing the necessary thermal 
energy during one scattering event. Instead, it must diffuse through the energy 
landscape with the help of many scattering events ~ (Vq/J) 2 >> 1. Since the cor- 
responding scattering matrix elements are suppressed as high powers of T /To and 
J /To, Ref. 19 , such diffusion must be strongly suppressed. The detailed solution 
of this problem along the line of the approach 23 will be presented elsewhere. 

The inconsistency of the boiling off mechanism has already been noted in 
Ref A Here we have shown that the alternative scenario is the fluctuative partial 
depinning of dislocation from the Peierls and 3 He pinning potentials. Within this 
scenario the corresponding activation energy seen experimentall y 24 i 25 i 15 should 
be attributed to the kink-pair creation energy (dependent on strengths of the both 
potentials) rather than to the activation of 3 He from a stationary dislocation. 

3 The puzzle of "missing dissipation" 

Both phenomena — the NCRI and the shear anomalies — are characterized by the 
dissipation peak in the region of fast variation of G(T) and the NCRI fraction 
Ps(T). The single-time relaxation scenario for such peak appears to be inade- 
quate 26 because the relaxation induced change of the real part appears to be signif- 
icantly larger than it should follow from the observed maximum of the imaginary 
part of the generalized response. This situation can be referred to as the puzzle 
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of "missing dissipation". In order to fix the problem on empirical level, various 
wide distributions of relaxation times can be invoked 2715 . As we have mentioned 
above, such conjectures, while providing a convenient empirical framework, do 
not explain the origin of the wide distribution in relatively good quality crystals. 

Here we argue that, in addition to the dynamical variations of the real and 
imaginary parts of the response which are connected by the Kramers-Kronig rela- 
tion, Peierls and 3 He pinning potentials induce a purely equilibrium temperature- 
variation of the real part (at zero frequency ft)) 14 which is not directly linked to any 
dissipative process. This relaxes the "constraint" between the real and imaginary 
parts of G. 

In order to illustrate the two origins of the shear modulus softening we intro- 
duce a linearized version of the model (QJ which takes into account the pinning gap 
A: all the non-linear terms are replaced by ~ K~ l A 2 y 2 /2. The corresponding dy- 
namical equation of the string biased by a small external stress a(t) ~ exp(— ift)/) 



with the boundary condition y (x — ±L /2 , t ) — 0, in addition to the standard relaxation- 
dissipation term ~ £ (see in Ref.—), includes the gap term ~ A 2 induced by Peierls 
and 3 He pinning. As shown above, A can be found from purely thermodynamical 
simulations of the full model ([T|) (see also in Ref.—). 

The nature of the dissipative kernel e(f') at low T represents a separate fun- 
damental problem. In this regard, we note the proposal 28 that dislocation kinks 
initiate strong dissipation in solid 4 He 29 . Similarly, the origin of the so called 
Bordoni peak in the internal friction in metals has also been assigned to kinks^i. 
Both approaches rely on the standard treatments in terms of the kinetic equation, 
without, however, addressing explicitly the issue of the integrability of the Sine- 
Gordon model (that is, that there should be no dissipation in it!). 

Here we conjecture that the main source of the dissipation are the kink-kink 
collision processes described in Ref. 31 under the assumption of the integrabil- 
ity broken by the spatial discretization of the Sine-Gordon description. Then, the 
long-time kink motion is described by diffusion coefficient D ~ exp(A/T)/A, 
RefM-. This corresponds to e(/') ~ 8(t') /D = C a A exp(-A /T) in Eq.© (where 
C ~ 1 in chosen units), so that the space-dispersive relaxation is given by ~ Dq 2 
with q standing for wave-vector along the dislocation. Such form implies that the 
dissipation vanishes as T — > 0, when the kink density ~ exp(— A /T)—>0, and also 
as the gap A (T) — > at large T, when kinks overlap and loose their meaning. 

Now we calculate the generalized shear modulus G(ft)) in the frequency do- 
main as G" 1 = CT, 1 + (Atid/L) j (y/o)dx (see the explanation above Eq.©). 
Then solving Eq.(|7| in Fourier with the zero-boundary condition we obtain 



y{x,t)-V 2 y(x,t)+ f dt'£(t')y(x,t-t')+A 2 y(x,t) = e(t), (7) 



G- 1 (r,ft)) = c;/ + 



Q 2 (T,co) 




tanh[r2(r,ft))L/2] 
Q{T,co)L/2 



(8) 
(9) 



The constant Aq = An^L 2 (which is independent of LI) is determined by the total 
softening effect. Specifically, at T = and ft) = 0, Q — and G 1 - G~/ ~ 
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Fig. 4 (Color online) The real (solid lines) and imag inary (dashed lines) parts of the variation of 
the generalized shear modulus G given in Eqs. d8|9l normalized by the total variation of G from 
T = to T = oo. The family of four curves corresponds to four values of the dissipation strength 
C shown in the plot. The other chosen parameters are: zli = 0.1, Aq = 0.05, L = 100, ft) = 
0.003, AqG c / = 0.4 (in the dimensionless units discussed in the main text). 



l/L 2 -s>0forL>> 1, and at ft) = and large T when A =0: G _1 (°o,0) — G^ 1 = 
Aq/12. 



A particular form of A(T) fitting well the experimental data was found to be 
A(T) = Ai(l — exp(— Aq/T)) 2 , where Aq depends on the strength of the Peierls 
and 3 He pinning potentials as well as on n and A\ ^0.1 is a constant. We have 
used these expressions in Eqs.( |8l9l > and have calculated the real and imaginary 
parts of the shear modulus, Fig. |4) As can be seen, in contrast to the simplistic 
single -relaxation-time model 26 ' 15 , the variations of the real and imaginary parts 
are not restricted by the condition that the maximum of the imaginary part is given 
by one-half of the total variation of the real part. Instead, while the real part is 
weakly dependent on the dissipation strength C a and is mostly controlled by the 
T-dependence of the gap A, the imaginary part is strongly dependent on the value 
of the friction coefficient C , and, in principle, can be zero. 



4 Kink-controlled core superfluidity 

Here we will give details of the proposal 1 1 linking the shear modulus anomaly to 
the superfiuid response. Its main ingredient is a strong interaction between the core 
superfluidity^*^ and density of kinks along the core. We will take into account 
this effect within the coarse grained mean field description relying on the bulk- 
averaged SF order parameter \ff as well as on the bulk-averaged density of kinks 
The interaction free energy between the two is taken in the minimal form 
F[ nt — g\n k | lff\ 2 > 0, with g\ > being some phenomenological coefficient. Thus, 
the full bulk free energy density takes the form 



F = a'{T - T )\\if\ 2 + ^\y\ 4 + gl n k \\if\ 2 + A(T) ■n k + Tn k ln(n k /(n Q e)), (10) 
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where d > 0,g 2 > are the standard Landau-expansion coefficients with To cor- 
responding to the mean field SF-transition temperature. The term ~ A describes 
energy of thermal activation of kinks and the last one accounts for the entropy of 
the kink gas in the limit «no, where no stands for the maximum bulk concen- 
tration of kinks determined by the size of the Sine-Gordon soliton « 1/ \f~A and 
the density of dislocations n^. Here kinks are treated as Boltzmann particles. 

We note that, formally, the model d 1 Ob apart from the kink-terms coincides 
with the approach 32 . However, in contrast to Ref. 32 advocating a true supersolidity 
in 4 He, we view the formation of the SF order parameter as being due to a network 
of dislocations with superfiuid core s 33 i 9 i 10 . In this analysis we don't consider that 
some dislocations exhibit superclimb and the anomalous compressibility 10 . The 
coarse grained description of such network is a separate interesting problem. 



4.1 The Mean Field phase diagram 

The equilibrium of the system can be determined from the minimization 8F /8\j/* = 
0md8F/8n k = 0: 



[a'{T-To)+gin k ]\tf + g2 \w\ 2 V = 0, n k =n exp^ ^ T ' J ■ (11) 

The first equation can be formally written as a relative value of the condensate 
density p s {T) = \ \(f\ 2 as 

p (rl = ftg) = i_l_ Sink 2 

p( ' p,(o) To g2Ps (oy 

where p s (0) = a'To/g 2 and we took into account that n^iT = 0) = 0. 

Eqs. dl lb describe the mean field phase diagram featuring lines of II and I order 
transitions with the tricritical point separating them. It is convenient to introduce 
the following dimensionless quantities 

Cl = ^Zll^(±)>0, C 2 = ^exp^)>0, (13, 
gin \TJ g\no \T J 

for T < To- Then, the II order (mean field) transition takes place along the line 
C\ = 1 , C2 > 1 . Along this line F has only one minimum iff = for C\ < 1 and 
<//■ ^ for Ci > 1, with n k 7^ in both phases. As C 2 < 1, F can have two coexisting 
minima. That is, the I order transition occurs for C2 < 1, with the tricritical point 
located at C\ =Ci = \. The actual equation for the I order transition line (when 
both minima have equal free energies) cannot be solved explicitly. Its numerical 
solution is represented by thick red line in Fig. [5] 

The boundaries (spinodals) within which F has two minima, that is, where 
metastable superfluidity can exist, can be found analytically. These are given by 
C 2 < l,C s (C 2 ) < Ci < 1 where 

Q(C 2 )=C 2 lnQ-Y 0<C 2 <1. (14) 



4 + «i|vf 



11 



C,1.4- 




Fig. 5 (Color online) The mean field phase diagram of the model i fTTT i in terms of the variables 
foi l. Solid red and dashed dark yellow lines represent the I and II order phase transitions, respec- 
tively. The large blue dot denotes the tricritical point separating the II and I order transitions. The 
domains above and below the transition lines marked as SF and as "normal" correspond to the SF 
and the normal states, respectively. The thin pink curve indicates the spinodal C s {Cz), Ea.JT4Ti. 
The metastable superfluidity can occur above this line and below the pink dashed straight line 
Ci = l,C 2 <l. 




Fig. 6 (Color online) Blue dots: Relative shear modulus, Eq.J2U. with T c = Q.2\K, from 4 . Red 
line - relative superfluid density obtained from Eq.<12t. Green dashed line shows zero for p s (T) 
at T = T C . 



4.2 Similarity between static superfluid and mechanical responses 

Here we will discuss the striking similarity between the SF response p s (T) and 
G(T) following from the model introduced above (cf. in Ref. 4 ). It occurs along 
the II order line due to the strong interaction between kinks and superfluidity given 
by the term ~ g\ in Eq.JTOll. Let's consider the system far from the tricritical point, 
C2 » 1, and close to the II order line, C\ ~ 1. The condition C\ = 1 determines 
the mean field transition temperature T c as 
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If 

gl»0 

a' T 

which constitutes the condition of strong kink-superfluid interaction, one finds 



» 1, (16) 



In (gi «o /a' 7b) 



« 7b. (17) 



At this point it is important to mention that the temperature T c is rather a renor- 
malized value of 7b than the actual transition of establishing the 3D coherence. In 
other words, it is the onset temperature where the NCRI should appear in a sense 
of the Shevchenko state 33 (cf. the vortex-fluid model—). 

The condition that the system is far from the I order transition C2 >> 1 in 
combination with T c « 7b, Eq.(fT6ll. gives 

" /( » 1. (18) 



gia% 



To what extent this condition is satisfied in solid 4 He remains to be seen. 

Let us demonstrate the similarity between the mechanical and superfluid re- 
sponses in the limit dl6l) . In the limit dl8l) T should be dropped from the definition 
of Ci in Eq. (ll3b for < T < T c , so that the solution for the superfluid fraction 
Ps(T) = |y/| 2 ,Eqs.CED, 

can be formally rewritten in the form 

p,(0) n k (T c ) 

where n k {T) obeys Eq. dl II) and T c is given in Eq. dl7b . 

The shear modulus thermal softening discussed above in Sec[2]is determined 
by the averaged density of thermal kinks n k . Indeed, Eq.© relates the G(T) soft- 
ening to the fluctuations of the string displacement (y 2 ) which, in its turn, is de- 
termined by the density of thermal kinks n k . Qualitatively, one finds (v 2 ) °= % at 
least at low T where kinks are still well defined entities. Thus, the shear modulus 
G given by Eq.© can be rewritten in terms of the coarse grained kink density as 

\+An k {T) 

where Go = G(T = 0) ~ G e i and the prefactor A can be taken as a constant of T 
for all practical purposes because it does not contain any exponential. 

We introduce the relative change of the shear modulus with respect to its value 
G(T C ) at the transition temperature (where p s is zero): 

Ml ' G(0)-G(T c )- { ' 

Substituting G from Eq.(l20ll, we find 

u(T) = I (l H^y\ « P^l (22) 

l+An k (T)V n k (T c )J ~ p,(0) (1Z) 
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where the last line follows from Eq.dT9l) and the factor 1/(1 +An(T)) » 1 because 
it controls the total relative variation of the modulus within only 5-15% — . Thus, 
within such accuracy, our model naturally features the close similarity between 
the relative variations of the superfiuid and mechanical responses. 

It should, however, be mentioned that the above simplified model cannot be ap- 
plied to temperatures significantly above the gap value A (T = 0), that is, when the 
simple activation dependence n k <x exp(— A /T) does not fit well (see in Ref. 14 ) the 
actual experimental data Ref. 4 . Nevertheless, Eq. dl2b does not actually rely on the 
simple activation form of the kink density and will remain the same regardless of 
a particular form of the kink free energy. Thus, the actual form of n k at elevated T 
can be extracted from the experiment 4 with the help of Eq. d20h . Then, substituting 
the so found n k (T) into Eq. dl2b and using g\/d as the only adjustable parameter 
for chosen Tq = IK, we find T c = 0.21^ and obtain the graph Fig|6] It shows the 
data from Ref. 4 , / = 2000Hz, plotted together with the relative superfiuid density 
p s {T) obtained from Eq. dl2b . The similarity between the two responses is obvious. 

Concluding this section, we note that the energy A sets the scale for both — 
G(T), below which it stiffens, and for the onset temperature T c of the SF response. 
As discussed in Sec l2.ll adding 3 He increases A and, therefore, increases the SF 
onset temperature. 

4.3 Similarity between the dissipative responses. 

The similarity between the SF and mechanical responses considered above is not 
limited to the static case. Obviously, a strong coupling ~ g\ must induce the same 
dissipation in both subsystems. As a result, excitations in the SF network will 
show damping similar to the one introduced into Eq.©. Such excitations lead to 
a dynamical mass (and the moment of inertia) redistribution inside the TO cell. 
Accordingly, the TO response should show the dissipation as long as SF currents 
percolate along a non-uniform network which, in general, allows linear coupling 
between SF phase and rotation. We demonstrate the similarity between the dissi- 
pation (in the shear modulus and in the SF excitations) in the mean field equations 
for small oscillations of the SF order parameter. Using the hydrodynamic approx- 
imation \j/ = y/ps"^ for the conjugate variables p, </> in F, Eq. dlOl) . and the canon- 
ical commutation relations, we find Hamiltonian equation for p from d 1 0b . For 
as classical particles we use linearized kinetic equation (cf. — ) in the single-time T 
relaxation approximation which ignores spatial dispersion: h k = — \8F j dn^. This 



where po denotes the coarse grained SF density; stands for the quasi-equilibrium 
density of kinks given by Eq. dl lb . The small variation Sn^ e9 ' of induced by 
the time dependent part p' of p = po + p' can be found from minimizing the free 
energy dlOb in the linear approximation as: 



gives 



P - PoV 2 (g2p+gin k ) =0, n k 




(23) 




(24) 
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where nx 9 ' stands for the static value of the kink density found from the second 



equation in (fTTT) . A substitution of <W e? ' into Eqs.d23l) and further elimination of 
in Fourier gives the dispersion relation for the SF excitations 



C0 A 



2 2 

-4<r 



c 



1(0T 
1 -iCOT 



0. 



2 ( e l) 

8i n o 

Tg2 



(25) 



denotes speed of 



where q stands for the wave- vector and Cq = Po(g2 - g\ n o^ /T) 
first sound. 

The value of T can be related to the relaxation in Eq.© as T _1 ~ £. Thus, at 
T — > 0, where T ~ exp(4 /T) — > °°, and at T — > °°, where T ~ 1 /A (T) — > °°, the 
dissipation in Eq. (l25b vanishes along with the dissipation of the shear modulus. 



5 Discussion and perspectives 

We have introduced the alternative to the traditionally accepted "boiling off model 
of pinning gliding dislocations by 3 He impurities: Quantum and thermal fluctua- 
tions of dislocation shape lead to decoupling of dislocations from the Peierls and 
3 He-pinning potentials at temperatures determined by the energy of kink-antikink 
pair. 

The analysis of the shape of the shear modulus vs temperature at significantly 
different 3 He concentrations, Ref. 4 , and its comparison with the results of Monte 
Carlo simulations strongly indicate that the relaxation time for establishing equi- 
librium between 3 He impurities bound to dislocations and those which are free in 
the bulk is very large on any experimental time scale. This conjecture is supported 
by the estimates of the values of the trapping potential for 3 He and is consistent 
with the general nature of the narrow-band impuriton 19 . The formation of the two 
subsystems of 3 He atoms — free in the bulk and bound to dislocations — is most 
likely to occur during 4 He crystal growth. We find this problem as very important 
for theoretical and experimental studies. 

It is shown that the shear modulus softening consists of two contributions: 
first, from thermally suppressed gap and, second, from dissipation. In general, 
these two effects are independent from each other so that the puzzle of "missing 
dissipation" can be resolved without envoking a wide distribution of relaxation 
times in relatively good quality crystals. 

The minimal model leading to the similarity between the TO NCRI response 
and the generalized shear modulus has been developed. It is shown that the strong 
suppression of the superfluidity along dislocation core by dislocation kinks locks 
in both responses. In light of the recent work- we note that the data which are well 
above the red line in Fig.l, Ref. 7 , does not exclude the genuine NCRI anomaly 
whose similarity to the mechanical response is determined by the mechanism de- 
scribed above. 

The dissipation similarly seen in the TO and the shear modulus measurements 
is proposed to be attributed to the bath of dislocation kinks living along the dis- 
location network. The nature of such dissipation (its frequency and temperature 
dependencies) is linked to the fundamental question of how dissipation emerges 
in a network of Id systems and, therefore, deserves close attention. Of a specific 
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interest is studying the dissipation in the system of superclimbing dislocations, 
where the spectrum is not sound-like anymore 10 . 

Another interesting aspect is how the SF fraction and its excitations can be 
detected in the UMass-Sandwich-type setup 8 . Currently, the detected flow is d.c. 
and is in the overcritical (strongly non-linear) regime and, therefore, is not suitable 
for answering these questions. Thus, it is important to devise and implement the 
a.c. current sandwich-type setup in order to study the linear response of the SF 
network directly. 

Acknowledgements We are grateful to John Beamish for providing the experimental data and 
useful discussions. We also thank Nikolay Prokof'ev and Boris Svistunov for many stimulat- 
ing discussions. This work was supported by the National Science Foundation under Grant 
No.PHY1005527, and by a grant of computer time from the CUNY HPCC under NSF Grants 
CNS-0855217 and CNS - 0958379. 



References 

1. E. Kim and M.H.W. Chan, Nature 427, 225 (2004); E. Kim and M.H.W. 
Chan,Science 305, 1941 (2004). 

2. A. J. Leggett, Phys. Rev. Lett. 25, 1543 (1970). 

3. E. P. Gross, Phys. Rev. 106, 161 (1957); E. P. Gross, Annals of Physics 4, 
57(1958); A. F. Andreev and I. M. Lifshitz, Sov. Phys. JETP 29, 1 107 (1969); 
D. J. Thouless Ann. Phys. 52, 403(1969) G. V. Chester, Phys. Rev. A 2, 256 
(1970). 

4. J. Day and J. Beamish, Nature 450, 853 (2007); J. Day,0. Syshchenko, and J. 
Beamish,Phys. Rev. B 79, 214524(2009). 

5. V.L. Tsymbalenko, Sov. Phys. JETP 60, 537(1984). 

6. J. Reppy, Phys. Rev. Lett. 104, 255301 (2010); X. Mi, E. Mueller, J. D. Reppy, 
larXiv: 1109.68181 

7. J. R. Beamish, A. D. Fefferman, A. Haziot, X. Rojas, and S. Balibar, Phys. 
Rev. B 85, 180501(R) (2012). 

8. M. W. Ray and R. B. Hallock, Phys. Rev. Lett. 100, 235301 (2008); Phys. 
Rev. B 79, 224302 (2009). 

9. M. Boninsegni, A. B. Kuklov, L. Pollet, N. V. Prokof'ev, B. V. Svistunov, and 
M. Troyer, Phys. Rev. Lett. 99, 035301 (2007). 

10. §. G. Soyler, A. B. Kuklov, L. Pollet, N. V. Prokof'ev, and B. V. Svistunov, 
Phys. Rev. Lett. 103, 175301 (2009). 

11. D. Aleinikava, E. Dedits, A. B. Kuklov, D. Schmeltzer, arXiv:0812.0983 . 

12. S. Balibar, Physics 3, 39 (2010): |http://physics.apsTor g/articles/v3/39; Nature 
464, 176 (2010). 

13. J. Friedel, Dislocations, Pergamon Press, 1964. 

14. D. Aleinikava, E. Dedits, A. B. Kuklov, and D. Schmeltzer, Europhys. Lett., 
89,46002(2010). 

15. O. Syshchenko, J. Day, and J. Beamish, Phys. Rev. Lett. 104, 195301 (2010). 

16. A. Granato and K. Lucke, J. Appl. Phys. 27, 583 (1956); ibid. 789(1956). 

17. P. Corboz, L. Pollet, N. V. Prokof'ev, and M. Troyer Phys. Rev. Lett. 101, 
155302 (2008). 



16 



18. C. Zhou, J.-J. Su, M. J. Graf, C. Reichhardt, A.V. Balatsky, I. J. Beyerlein, 
arXiv: 1110.0841. 

19. A.F. Andreev, Sov. Phys. Usp. 19, 137(1976). 

20. N.V. Prokof'ev and B.V. Svistun ov, Phys. Rev. Let t. 87, 160601 (2001). 

21. D. Aleinikava and A.B. Kuklov. larXiv:1108.1182l 

22. A good fit can only be achieved for the very unrealistic values: ~ 
10- 2 ,X 3 - W~ 2 ,E a = 0.2-0.3K. 

23. Yu. Kagan and LA. Maksimov, Sov. Phys. JETP 38, 307(1974). 

24. I. Iwasa, H. Suzuki, J. Phys. Soc. Japan, 49, 1722(1980). 

25. M. A. Paalanen, D. J. Bishop, and H. W. Dail, Phys. Rev. Lett. 46, 664 (1981). 

26. C.-D. Yoo and A. T. Dorsey, Phys. Rev. B 79, 100504(R) (2009). 

27. Z. Nussinov, A. V. Balatsky, M. J. Graf, and S. A. Trugman, Phys. Rev. B 76, 
014530 (2007). 

28. A.V. Markelov, Sov. Phys. JETP 61, 118(1985). 

29. Y. Hiki and E Tsuruoka, Phys. Lett. A 62, 50(1977); F. Tsuruoka and Y. Hiki, 
Phys. Rev. B 20, 2702(1979); V. L. Tsymbalenko, Sov. Phys. JETP 47, 1507 
(1978). 

30. F. Marchesoni, Phys. Rev. Lett. 74, 2973(1995). 

31. K. Damle and S. Sachdev, Phys.Rev.Lett. 95, 187201 (2005). 

32. PW. Anderson, Science 324, 631 (2009); 

33. S. I. Shevchenko, Sov. J. Low Temp. Phys. 13, 61 (1987); D. V. Fil and S. I. 
Shevchenko, Phys. Rev. B 80, 100501 (2009). 

34. P. W. Anderson, Nature Physics 3, 160 (2007) 



